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1 Introduction 

Nowadays there are several experimental results that may hint to a deviation from the 
coherent picture of the standard three neutrino oscillation model. These results, not stati¬ 
stically signihcant on a single basis, could be explained with the hypothesis of the existence 
of additional neutrino mass eigenstates at eV mass scale, namely sterile neutrinos [T]. In 
particular there are strong tensions between Vf, disappearance and Ve appearance results 
in experiment at Short-Baseline (SBL). The main source of tension arises from the ab¬ 
sence of any disappearance signal [2]. Limited experimental data are available on 
disappearance at SBL: the dated CDHS experiment [3] and the more recent results from 
MiniBooNE [1], a joint MiniBooNE/SciBooNE analysis [5] and MINOS and SK [6]. The 
latter results only slightly extend the disappearance exclusion region set by CDHS and cur¬ 
rently the region s\v?29new < 0.1 is largely still uncostrained. 

After the decision of CERN not to pursue for the time being the construction of a new 
neutrino beam, the NESSiE Collaboration proposed to perform an accurate measurement 
of disappearance at small L/E using magnetic spectrometers (B = 1.5 T) at NEAR 
and EAR sites on the ENAL-Booster neutrino beam [an]. Many detector configurations 
were considered investigating experimental aspects such as the measurements of the lepton 
charge on event-by-event basis and its energy over a wide range. The basic NESSiE con¬ 
cept is to design, construct and install two spectrometers at two sites, “Near” (at 110 m, 
on-axis) and “Ear” (at 710 m, on surface), in line with the ENAL-Booster, fully compatible 
with the proposed LAr detectors. Using very massive detectors (300 tons and 700 tons 
for the Near and Ear spectrometers, respectively) it would be possible to collect very large 
statistical samples in order to improve the disentangling of systematic effects and to span 
the oscillation mixing parameter down to till un-explored regions {s\.v?29new > 0.01 or less). 
By keeping the systematic error at the level of 1-2% for the measurements of the inte¬ 
ractions, it would be possible to: 

• measure disappearance in the almost entire available momentum range (p^ > 500 
MeV/c). 

• measure the neutrino flux at the Near detector, in the relevant muon momentum 
range, to keep the systematic errors at the lowest possible values; 
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• measure the sign of the muon charge to separate from for systematics control. 

The uncertainty on the absolute flux at MiniBooNE stays below 20% for energies below 
1.5 GeV while it increases drastically above that energy. The uncertainty is dominated 
by the knowledge of hadronic interactions of protons on the Be target, which affects the 
angular and momentum spectra of neutrino parents emerging from the target. The result 
is based on experimental data obtained by the HARP and E910 collaborations [9]. 

Such a large uncertainty makes the use of two or more identical detectors at different 
baselines mandatory when searching for small disappearance phenomena. The ratio of 
the event rates at the Ear and Near detectors (FNR) as function of neutrino energy is 
a convenient variable since it benefits at first order from cancellation of common proton- 
target and neutrino cross-sections systematics and of the effects of reconstruction efficiencies. 
Thanks to these cancellations the uncertainty on the FNR or, equivalently, on the Far 
spectrum extrapolated from the Near spectrum is usually at the percent level ranging in 
the 0.5-5% interval. In order to understand how the hadroproduction uncertainty affects 
the accuracy on FNR for the specific case of our proposal we developed a new beamline 
simulation using either FLUKA, GEANT4 or the Sanford-Wang parametrization for the 
simulation of p-Be interactions [9] . 

2 Physics analysis and sensitivity 

We developed sophisticated analyses to determine the sensitivity region that can be explored 
with an exposure of 6.6x10^^ p.o.t., corresponding to 3 years of data taking on the FNAL- 
Booster beam. Our guidelines were the maximal extension at small values of the mixing 
angle parameter, and its dependence on systematic effects. To this aim the sensitivity of the 
experiment was evaluated using three different analyses implementing different techniques 
and approximations: 

• a Feldman&Gousins technique with ad hoc systematic errors added to the muon mo¬ 
mentum distribution; 

• a full correlation matrix based on the full Monte Carlo simulation including the re¬ 
construction of the simulated data; 

• CLs profile likelihood; 

For all analyses the two flavor neutrino mixing in the approximation of one mass dominance 
is considered. The oscillation probability is given by the formula: 

P = sin^ {26 new) sin^ {1-27 (km)/E {GeV)) (1) 

where Arrinew is the mass splitting between a new heavy-neutrino mass-state and the heav¬ 
iest among the three SM neutrinos, and Onew is the corresponding mixing angle. 

In the second method we implemented different smearing matrices for two different obser¬ 
vables evaluated with GLoBES m, the muon range and the number of crossed planes at 
Near and Far detectors respectively, associated with the true incoming neutrino energy. 
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These matrices were obtained through a full Monte Carlo simulation. In this analysis the 
correlations between the data collected in the Far and Near detectors are taken into account 
through the covariant matrix of the observable. 

The disappearance can be observed either by a deficit of events (normalization) or, alter¬ 
natively, by a distortion of the observable spectrum (shape) which are affected by systematic 
uncertainties expressed by the normalization errors matrix and the shape errors matrix, re¬ 
spectively. The normalization errors matrix is the component of error matrix which is the 
same for each element. The shape errors matrix represents a migration of events across the 
bins. 

By applying the frequentist method the statistic distribution was used to determine the 
sensitivity to oscillation parameters. Sensitivity plots were computed by introducing bin-to- 
bin correlated systematic uncertainties as expressed in the covariance matrix by considering 
either 1% correlated error in the normalization or 1% correlated error in the spectrum shape. 
As a representative result the sensitivity calculated considering 1% correlated error in the 
shape is plotted in Fig. [TJ The sensitivity computed considering CC and NC events is 
almost the same as the sensitivity obtained with only CC events m and therefore NC 
background events do not affect the result. 



Figure 1: 95% C.L. sensitivity obtained using the muon range for all CC (black) and 
CC-I-NC (red) events and for only CCQE events (blue) assuming 3 years of data taking. 
1% bin-to-bin correlated error in the shape is considered. 


3 Summary 

Existing anomalies in the neutrino sector may hint to the existence of one or more addi¬ 
tional sterile neutrino families. The NESSiE collaboration performed a detailed study of 
the physics case in order to set a Short-Baseline experiment at the FNAL-Booster neutrino 
beam to exploit the measurement of the charged current events. 
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The best option in terms of physics reach and funding constraints is provided by two spec¬ 
trometers based on dipoles iron magnets, at the Near and Far sites, located at 110 and 710 
m from the FNAL Booster neutrino source, respectively. 

By the proposed approach systematic errors can be kept at the level of 1-2% for the mea¬ 
surements of the interactions. In particular the measurement of the muon momentum 
can be performed at the percent level and the identification of its charge on event-by-event 
basis can be extended to well below 1 GeV. 
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